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Abstract 



We present results of numerical simulations of the formation of black holes from the gravitational 
collapse of a massless, minimally-coupled scalar field in 2+1 dimensional, axially-symmetric, anti de-Sitter 
(AdS) spacetime. The geometry exterior to the event horizon approaches the BTZ solution, showing no 
evidence of scalar 'hair'. To study the interior structure we implement a variant of black-hole excision, 
which we call singularity excision. We find that interior to the event horizon a strong, spacelike curvature 
singularity develops. We study the critical behavior at the threshold of black hole formation, and find 
a continuously self-similar solution and corresponding mass-scaling exponent of approximately 1.2. The 
critical solution is universal to within a phase that is related to the angle deficit of the spacetime. 

1 Introduction 

The past several years has seen growing interest in the properties and dynamics of asymptotically anti de- 
Sitter (AdS) spacetimes, predominantly due to the discovery of black hole solutions in 2+1 dimensional AdS 
spacetime 0] and the AdS/CFT conjecture The existence of vacuum [] black holes (also called BTZ 
black holes) is surprising because the local solution to the field equations is isometric to AdS, and hence 
has constant curvature. What makes a BTZ spacetime different from AdS is its global structure, which 
can be obtained by making appropriate identifications within the universal covering space of AdS Q. The 
natural question that such a construction poses is: how similar are these black holes to their more familiar 
3+1 dimensional (4D) counter-parts? In particular, do these black holes have thermodynamic properties 
when considered within the framework of a quantum theory, and can they form through dynamical collapse 
processes? It turns out that BTZ black holes do bear striking resemblance to 4D black holes in many respects 
(see Q for review articles). In this paper we present the results of a numerical study of the collapse and 
formation of non-rotating BTZ black holes from a massless scalar field in 2+lD AdS spacetime. Of particular 
interest is whether critical phenomena |^ are present at the threshold of black hole formation — namely if 
by fine-tuning of initial data, we can make the system asymptote (at "intermediate times") to a solution 
which is universal in the sense of not depending on details of the initial data. Furthermore, if the black hole 
transition is "Type 11" , so that there is no smallest mass of black hole which can be formed, then we expect 
their to be a scaling relation for the black hole mass of the form M = K{p — p*)^^ . Here p is a parameter 
in a family of initial data such that p — p* is the critical solution, X is a family dependent constant and 7 
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is a universal exponent (see p| for a recent review). The 'extra' factor of 2 in the exponent is expected for 



BTZ black holes — see section 3.2. As we will show, it turns out that the system does exhibit a continuously 
self-similar (CSS) solution in the critical limit, with a scaling exponent 7 = 1.2 ± 0.05. 

Earlier works on black hole formation in AdS considered disks of dust Q, null radiation thin dust 
rings 1^ , and the collision of point-particles |0 . In the case of dust- ring collapse Peleg and Steif found a 
scaling exponent of 1/2 at the transition between black hole and naked singularity formation. Birmingham 
and Sen found the same exponent at the threshold of formation in the case of colliding particles. Husain 
and Olivier have also studied the massless scalar field in 2-t-l dimensions using a double null formalism, and 



have formed black holes with their code |11 



Our paper is organized as follows. In section |^ we describe the system of coordinates and numerical 
scheme we have chosen to use, and the resultant field equations and boundary conditions. An interesting 
consequence of our analysis is that we are unable to derive boundary conditions for the scalar field at the 
edge of the universe that are analogous to the out-going radiation conditions often employed in numerical 
relativity. In AdS spacetime the scalar field reaches time-like infinity 2 in finite proper time as measured 
by a central observer, and the only consistent boundary conditions we can place on the scalar field confine 
it to the universe. This is reassuring from the standpoint of global energy conservation, but complicates 
the search for the universal scaling relation between black hole mass and parameter-space distance to the 
critical solution. The system behaves as if the scalar field is within a finite sized box, and so when a black 
hole forms all of the scalar field initially present eventually falls into the hole. M{p) is therefore trivially a 
function of how the initial energy distribution scales with p. 

In section ^ we present results from the evolution of several families of initial data, focusing on critical 
behavior. To obtain 7, we follow the work of Garfinkle and Duncan [p^, and examine the scaling of the 
maximum value attained by the curvature scalar R in the sub-critical regime. We also study the effect that 
a central point particle (characterized by the angle deficit of the spacetime) has on the critical solution. As 
expected, we find that the more massive the point particle, the smaller the initial amplitude of the scalar 
field that gives rise to the critical solution. One might thus expect to have a one-parameter family of critical 
solutions with an overall scale related to the particle mass. It is surprising, therefore, that the scalar field 
always grows to the same amplitude in a near-critical evolution. A phase shift in central proper time is the 
only qualitative difference attributable to the mass of the particle. At the end of section || we study the 
interior structure of black holes that form, giving evidence that a 'crushing' spacelike curvature singularity 
forms within the event horizon. Thus the interior structure is significantly different from the BTZ solution, 
which has constant curvature (though the BTZ singularity is still crushing for extended objects falling into 
it). 



2 The Einstein Klein-Gordon system in AdS spacetime 

We solve the Einstein field equations in 3 spacetime dimensions with cosmological constant A = — 1/^^, 
coupled to a massless Klein-Gordon (KG) field 

Rab - -^Rgab + ^gab = uTab, (1) 

where the stress-energy-momentum tensor for the KG field cj) is [|l^ 

Tab = (t>;a4lib - -^gab4';c4'''' ■ (2) 

Covariant differentiation is denoted by a semi-colon, while a comma denotes partial differentiation. We 
only consider circularly symmetric configurations of a minimally-coupled scalar field in this paper. Hence, (j) 
satisfies the wave equation 

□0 = = 0, (3) 

and in coordinates {t, r, 6) adapted to the symmetry, characterized by the Killing vector d/d9, (j){r, t) is only 
a function of the radial coordinate, r, and time coordinate, t. 
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One of the many peculiar features of AdS spacetime is its causal structure. In particular, null infinity 
I is time-like, and any observer living in AdS spacetime can send and receive light-like signals to and from 
X in finite proper time p4|. These properties of AdS make it challenging to deal with numerically, as the 



scalar field traverses the entire universe on a local dynamical time-scale. Also, as we will show in section 2.1 
the only regular boundary conditions on the field (/> at X are Dirichlet conditions, so we cannot ignore the 
unusual causal structure of the spacetime by, for instance, placing out-going radiation boundary conditions 
on at a finite proper distance from the origin. For these reasons, we adopt a coordinate system in which 
the metric takes the form: 

t) and B[r,t) are arbitrary functions of (r, t), and it is straight-forward to show that when A = B = 
the above metric describes AdS spacetime; i.e. it is a solution to ([l]) with Tab = 0. Notice that, in this metric, 
radial null geodesies travel with constant coordinate speed dr/dt = ±1, and X is at r = 7r£/2. The metric 
is singular at X, but we can place regular boundary conditions on A and B there, so that the spacetime 
is asymptotically AdS. Also, if we interpret as a periodic angular variable then the above metric has the 
correct topology to represent a BTZ black hole, as the topological censorship theorems require that the 
boundary at infinity share the topology of any event horizon that may exist in the interior of the spacetime 
[ p^ . However, for the non-rotating collapse described in this paper, 9 has no dynamical significance. 
Defining 

<f(r,t) = 0,, , n(r,0-<^,t (5) 
and using units where k = 47r, we get the following set of equations upon expanding (|l|)-(^ with the metric 

A..-A« + ^l^^ + 2.(a>^-tf).o, (6) 

B,rr - B,, + B,. i^B,r + j—^-j^^—^^ - {B,f + = 0, (7) 

B.. + B.IB.-A^ ■ l + cos>/.) 



icos{r/£) sin(r/^) 



.cos(./^);in(./f) - + Zfe^ + '""^^^ + = 

+ (^'^ - + ^) - + .sin(./.?cos(./J + ^-'^n - ' 

and 

[tan(r/^)e^$] ^ - tan(r/^) [e^H]^ = 0. (10) 

Within the context of the 3+1, or ADM, formalism, equations (§) and are the Hamiltonian and 
momentum constraints respectively, while equations (^) and are combinations of the evolution and 
constraint equations. Equation ( p^ ) is the wave equation for the scalar field. There are two unknown 
geometric variables — A{r, t) and B{r, t); hence one needs to use at least two of the four equations ^ - (^) to 
dynamically determine the geometry. In this work, we have chosen to use equations (1) and (0) to update A 
and B. As is common practice in such a "free evolution scheme" , we can then use residuals of the constraints 
(H) and as one way of estimating the level of error in our solution. 

With regards to initial conditions, we choose to freely specify $(r, 0) and n(r, 0) (we have to specify two 
scalar-field degrees of freedom at each r), as well as B{r, 0) and i3,t(r, 0). A{r, 0) and A_t(r, 0) are then fixed 
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from the constraint equations (see sec. 2.2 for more details). This procedure is clearly somewhat ad hoc, but 
has worked very well in our study. 

The Ricci scalar of this spacetime is 



R 



($2 



6 



(11) 



The Weyl tensor is zero, and other non-zero curvature scalars can be expressed as polynomial functions of 
R. 



2.1 Regularity conditions 

We require that the solution for our dynamical variables A{r,t), B{r,t),^{r,t) and Il{r,t) be regular at 
the origin, r — 0, and at X, r = tt£/2. The field equations then essentially dictate the allowed boundary 
conditions on these variables. By inspection of (|^)-([lO|) we obtain the following conditions. At r = 



and at r = tt£/2 



A{TT£/2,t) 

B^r{TT£/2,t) 

$(7r£/2,<) 
U{7r£/2,t) 



AAO.t) 

A^r{0,t) 

B,(0,i) 
$(0,t) 

n.(o,i) 
= A. 



0. 



St(0,i) 





{TT£/2,t)=A,tiTT£/2,t)^0 



(12) 
(13) 
(14) 
(15) 
(16) 

(17) 
(18) 
(19) 
(20) 



Note that condition (|l^ ) on 11(0, t) is a direct consequence of the defining relation for Il{r,t) and the 
regularity condition for $(0, t) (|l5|). Also note that we have multiple conditions for B at the outer boundary, 
and for A and B at the origin. We have chosen to implement the Neumann conditions for A and B at the 
origin and the Dirichlet condition for A at X, and then to monitor the other conditions as a consistency 
check during evolution. Conditions (pT|)-(|2C|) ensure that the spacetime is asymptotically AdS. 

It is interesting that the field equations enforce Dirichlet boundary conditions on $ and H, effectively 
preventing us from implementing out-going radiation boundary conditions at X (if we wanted to let the field 
"leak out of the universe" when it reaches T). To see this more clearly, consider the energy fluxes Tabri'^'r]^ 
and Tab£°'£^ along outgoing and ingoing null vectors, £"" and 77", respectively, normalized so that £°'ria = — 1 



t = 



cos{r/£) 
~J2^ 



d_ 

dt 



d_ 

dr 



(21) 



_ cos(r/^) 



^/2e^ 



dt dr 



A straight forward calculation using (0) gives 



E± = 



cos(r/^)2($±n)2 



2e 



2A 



(22) 



(23) 



where E+ is the influx and E- the outflux. Thus no- outflux/influx boundary conditions can be obtained 
in the usual way by differentiating $ ± 11 with respect to r and t in turn, and utilizing the fact that, from 

(|, a>,t = n,,: 

$,^±$,4 = (24) 

n,^±nt==o. (25) 

(26) 
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Here, the plus sign corresponds to no-influx, and the minus sign to no-outflux. However, at the outer 
boundary, regularity forces $(7r£/2,i) = U{7r£/2,t) = 0, and hence $,t(7r£/2,t) = n_t(7r€/2,i) = 0, so there 
is no distinction between the no-influx and no-outflux condition. The only situation consistent with both 
conditions is that no flux crosses the outer boundary in either direction. Even when we try to derive no- 
outflux/influx conditions with the asymptotic behavior of (j) factored out, namely letting = cos^ {r / 
and placing boundary conditions on ^, we find that the wave equation on T cannot distinguish between 
no-outflux and no-influx conditions. Also, in early experiments we were unable to obtain stable numerical 



evolution with the no-influx boundary conditions (|24| ) and (25) applied at a flnite proper circumference, 
corresponding to r < t:1/2. The Dirichlet boundary condition at T is also consistent with the behavior of 
a massive scalar field in an AdS background, where an infinite effective-potential barrier prevents any of 
the field from reaching X, regardless of how small the mass is. Of course, all of this does not mean that an 
effective outgoing radiation condition can not be implemented for the massless field in asympotically AdS 
spacetimes. In any case, in the context of the current study, we would be apt to view such a condition as a 
numerical convenience, rather than being of any intrinsic physical interest. 

2.2 Initial conditions 

For initial conditions at t = 0, we are free to specify the scalar field gradients <i>(r, 0) and n(r, 0), the metric 
function B{r, 0) and its time derivative B t(r^ 0). We then numerically solve for A{r, 0) and A.t(r, 0) using the 
hamiltonian and momentum constraints and (^. The freedom that we have to specify i?(r, 0) amounts 
to a choice of what the proper circumference, (€tan(r/£)e^), and its initial time derivative are, as a function 
of the radial coordinate r. That we do not have the freedom to choose B for all time is a consequence of 
the gauge condition that radial light-like signals travel with unit coordinate velocity. For simplicity we set 
B(r,0) = B^f{r,0) = 0. 

We believe (though are unable to prove so), that the set of conditions just described is capable of 
generating all possible initial data, which is regular and free of trapped surfaces, for the minimally-coupled 
scalar field in asymptotically AdS spacetime (in 2-1-1 dimensions). The presence of trapped surfaces at t = 
is incompatible with the conditions on B{r,0) and B,t{r,0) — in our coordinate system dr/dt = 1 along 
an outgoing null curve, and hence a non-zero _B(r, 0) and/or i?_((r, 0) is required to describe non-positive 
outward null-expansion. However, in this study we are only interested in initial data that is free of trapped 
surfaces, so the conditions on _B(r, 0) and B^t{r,Qi) are not restrictive. 

For the initial scalar field profile, 0(r, 0), we consider three families of functions — a gaussian curve raised 
to the n*'' power 

(/)(r,0) =Pe(('-''°)/'^)'", (27) 
a 'kink' (based on an arctan function) for which $ = d(j)/dr is 

-2PVo'cos(r/^) sin(r/£) [-^ sin(r/^) cos(r/£) 2(r - ro)(l - 2 sin(r/€)2)] e-('-'~")'/'"' 

[a sin(r/^)4 cos(r/^)4 -f (r - tqY] 

and a family of harmonic functions ^ 

0(r,O) = Pcos2(rn/£), (29) 

where P,rQ,a and n are constant parameters. Then, depending upon whether we want to model initially 
ingoing, outgoing or static fields, we set n(r, 0) = $(r, 0), n(r, 0) = — <i>(r, 0) or n(r, 0) = respectively. 
Note that this method cannot give purely ingoing or outgoing pulses — n(r, t) = ±$(r, t) is not an exact 
solution to the wave equation, and a little bit of energy always propagates in the opposite direction to that 
desired. 

As noted previously, we set B{r,0) — _B j(r, 0) = 0. The remaining geometric variables, A{r,0) and 
A t{r,0) are then computed from the Hamiltonian and momentum constraints (||) and (|9|). We integrate the 
constraints outwards from r = 0, setting A^t(0, 0) = 0. For the most part we will consider the collapse of 

■^we call these functions 'harmonic' because without back-reaction and for initially static configurations (n(r, 0) = 0) the 
exact solution to the wave equation is periodic in time. 
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a scalar field initially exterior to empty AdS space. This corresponds to setting A{0, 0) = 0. However, in 
section 3.2.2, we will briefly consider the effect of collapsing the field in the presence of a point particle at 



the origin, the calculation of which involves introducing an angle deficit into the spacetime. From the metric 
(or by examining the parallel transport of a vector about r = in an infinitesimal loop), the angle deficit lu 
at i = is related to A{0, 0) as follows: 

c^ = 27r(l-e^(''^°)). (30) 

Of more interest is the relationship between A{0, 0) and the mass of the point-particle, Mpp. The remainder 
of this section is devoted to finding this relationship, and in the process we will define a general mass aspect 
function M{r,t) for the spacetime. 

When the scalar field gradients identically vanish (which they do at T, and, to an excellent approximation, 
at r = for the initial data that we consider) , the Hamiltonian constraint has the simple solution 

k-cos^{r/e)' ^ ' 

where fc is a constant of integration. We can relate k to the BTZ mass parameter M of the spacetime by 
appealing to the usual form in which the BTZ solution is expressed: 

ds^ = -{-M + f^/e)dP + —TT^^-^^df'' + f'^de'^. (32) 

— Ai + r'^ 1 

M = —1 is AdS spacetime, M > are black hole solutions and M < 0, M ^ — f are spacetimes with conical 
singularities, or point particles at the origin (the range of f is from to oo). For general (non- vacuum) 
solutions let us define the mass aspect M(r,t) as follows 

I Vrp EE -M(f, t) + f^jl^. (33) 

Then, in our coordinate system (^, M{r,t) takes the following form 

M{r, t) = e2(^-^) [e^^ tan2(r/^) + f sn?[r/(.)[[B^t? - (S r)') - 2£tan(r/£)B,, - se^{r/i)\ . (34) 

Using the field equations (^)-(||) it is straightforward to show that M is a conserved quantity in regions of 
the spacetime where $ and 11 are zero (in particular at X). At t — 0, where B = B^t — 0, we can substitute 
"n]) into (H) to find k: 

^ (35) 



1 + M 



When M > Q our metric ^ with the chosen initial conditions is singular at the horizon of an empty BTZ 
spacetime, but the metric turns out to be well behaved at < = for initial data that does not contain trapped 
surfaces 0. 

Finally, from (|3l|) and (|3^) the contribution, Mpp, of the point particle to the mass of the spacetime is 

Mpp = f - e-^^^^'O) (36) 



2.3 Numerical Scheme 

We solve the set of equations (|^),(0) and (^0|) by converting them to a system of finite difference equations 
on a uniform coordinate grid using a two-time level Crank-Nicholson scheme. We also add Kreiss-Oliger 
style dissipation p6| to control high-frequency solution components; this is crucial for the stability of our 
method. 

At first, we used standard 2nd order accurate 3-point finite difference stencils for the spatial derivatives 
at each time level — centered-difference operators at interior points, a forward-difference operator at the inner 

''however, because of our choice of gauge, we know that the coordinate system must become singular within one hght-crossing 
time (LCT) of the formation of a black hole. The event horizon is a null hypersurface travelling outward with unit coordinate 
velocity, thus the coordinate distance between the event horizon and I will go to zero within a time t = ■Ki/2. 
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boundary and a backward-difference operator at the outer boundary. However, we found that these operators 
excited a small instability in the metric variables in the vicinity of the outer boundary. The resultant ripples 
would propagate inwards and cause problems in situations where black hole formation was imminent. The 
primary source of these ripples was truncation error in the solution A exciting small oscillations in B. 
Specifically, A acts as a source term in the evolution equation for B (|^) , and B happens to be very sensitive 
to small errors in A near the outer boundary (essentially since the leading term of A, when considered as a 
power series in cos^(r/£), cancels with the spatial derivatives of B in (|^) initially, and so higher order, less 
accurately known terms of A are responsible for i?'s "acceleration"). To reduce these problems we now use a 
5-point, 4th order accurate spatial derivative operator at interior grid points, and 6-point 4th order backward 
and forward operators near boundaries that have the same truncation error as the interior operator. Also, 
we find that using the momentum constraint to solve for A at the next to last grid point is necessary to 
obtain convergence of the solution as we go to finer spatial resolution (for some as yet unknown reason the 
evolution equation was exciting a growing mode on finer grids at that point). The program to perform the 
evolution was written in Fortran 77 and RNPL (Rapid Numerical Prototyping Language ||l^); animations 
and pictures from several evolutions can be obtained from our website jl^. 



2.4 Detecting black holes and excising singularities 

To detect black hole formation we search for trapped surfaces, defined to be surfaces where the expansion of 
outgoing null curves normal to the surface is negative. If cosmic censorship holds, then trapped surfaces are 
always found within the event horizon of a black hole, though at the end of the simulation we can trace null 
rays backwards from 2 to confirm this. In our coordinate system the condition for a surface to be trapped is 

1 + ecos{r/e) sin{r/£){B^r + Bj) < 0. (37) 

We estimate the mass of the black hole by monitoring the proper circumference 27r^ tan(rAH/^)e^^''*'^'*-' of 
the apparent horizon (the outer-most trapped surface), and use the relationship between BTZ black hole 
mass and event horizon circumference ((^) — the horizon is at f = y/M£)): 

M « tan2(rAH/^)e2^^"^'''*^ (38) 

If all of the scalar field is absorbed by the black hole during evolution, then the estimated mass should 
eventually become equal to the initial, asymptotic mass of the spacetime as given by ( p^ ) in the limit 
r ne/2. 



As we will show in section 3.3, shortly after an apparent horizon(AH) forms, we find what appears to 
be a spacelike curvature singularity forming within the AH. If we use a straightforward evolution scheme, 
the metric and scalar field variables quickly diverge, and any given simulation just as quickly breaks down. 
At the same time, we would like to probe the structure of the spacetime approaching the singularity, as 
well as to continue to following the evolution outside the AH as long as our coordinate system allows 
(approximately 1 light-crossing time). To accomplish this, we have implemented singularity excision, a 
technique fundamentally motivated by the black- hole-excision strategy first proposed by Unruh [|l^ . 

Our excision strategy is as follows. We monitor the magnitude of the metric variables, and when they 
grow beyond a certain threshold Q at any point we excise that point plus a small buffer zone (of 4 to 6 grid 
points) on either side of it. (Note that the non-excised region of the grid will no longer be contiguous if the 
excised point is further away from the original grid boundaries than the size of the buffer zone) . At the new 
grid boundaries exterior to the excised region, we continue to solve for the metric and field variables using 
the evolution equations, but replace all centered-difference operators with forward and backward-difference 
operators, as appropriate, so that the solution is not "numerically influenced" by the excised grid points. 
Physically, the solution that one would obtain within the causal future of the excised zone is meaningless, 
so we also remove this region of the grid during subsequent evolution. In our coordinate system this is easy 
to implement, as radial null curves travel at constant, unit coordinate velocity. Thus, if our grid-spacing is 
Ar, after an amount of time At = Ar we expand the excised region by 1 grid point on either side. Also, 
we continue to monitor the metric variables on the remainder of the grid, and when they grow beyond the 

■^for a threshold we choose a number that is sufficiently large so we are fairly certain (from past experiments) that if any 
variable grows beyond the threshold then a crash is imminent 
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threshold at any other points we expand the excised region to include those points (and a buffer) . Thus the 
excised piece of the grid is always contiguous. In principle, it would not be difficult to keep track of multiple 
excised zones, though we did not find it necessary to do so for the interior solution shown in sec. I^a 
single zone is sufficient to obtain a good view of all of the interior up to the putative spacetime singularity. 

We have tested the singularity excision scheme by excising a light cone from a solution that remains 
regular, and verifying that the excised solution does converge to the regular solution as Ar decreases. 

In summary, we briefly clarify the difference between singularity and black hole excision. First, notice that 
we never use trapped surfaces to trigger the excision of a region of the grid. Thus, our code could, without 
modification, excise naked and coordinate singularities. The boundary of the excised region is always null or 
spacelike, so the scheme might not be able to distinguish between timelike and null curvature singularities. 
However, if a timelike singularity was encountered, it may still be possible to deduce its nature by examining 
the curvature invariants just exterior to the excised surface. For example, suppose during evolution a light- 
like region was excised, and curvature invariants started diverging as one approached the initial excised point, 
yet remained relatively "small" and finite just outside the future light-cone of the excised point, then one 
would have reasonable evidence for a timelike singularity. Second, with the singularity excision scheme, we 
excise only the region of the grid to the causal future of the singularity. In the case of a black hole spacetime, 
this results in a more complete view of the spacetime than what one would obtain with the standard black 
hole excision strategy (which would have in Fig. 19, for example, excised the region of the spacetime labeled 
"region of trapped surfaces", and everything to the left of it). 



3 Results 

In this section we discuss results from the evolution of several sets of initial data, focusing on the threshold 
of black hole formation. For convenience we set £ = 2/7r so that X is at r = 1, though the results presented 
here are valid for any non-zero, finite £, through an appropriate rescaling of the metric variables and scalar 
field gradients. Specifically, consider the following coordinate transformation 

with f defined on the range [0, 7r/2]. Then it is easy to see that the £ dependence cancels from all equations 
(p[)-(p^ when expressed in terms of r and i. So, given a solution A{r,t), B{r,i), $(f,i) and n(f, i) to the 
rescaled field equations we can find a corresponding solution for any £ by inverting the transformation ( p9| ) 
(see also (||)): 

A{f,i) A{r/£,t/£) (40) 

BifJ) ^ Bir/£,t/£) 
$(f,?) £^{r/£,t/£) 
n(f,i) ^ £Ii{r/£,t/£), 

with r ranging from to n£/2. Notice that the initial energy density, being proportional to ($^ + 11^), scales 
like so there is no straight-forward method to extrapolate a solution to the limit of zero cosmological 
constant, where £ ^ oo. 

We present results from 4 families of initial data: an ingoing gaussian ((p^) with n = 1), an ingoing 
squared gaussian ((p7|) with n = 2), an ingoing kink (^), and a time-symmetric, n — 1, harmonic function 
(P9|) . In each case we vary the amplitude P when tuning to the black hole threshold^, and for the first three 



families we have chosen a — 0.05 and tq = 0.2. Except in section 3.2.2, where we briefly study collapse onto 
a point-particle, we have set A{0,0) = in all cases, corresponding to angle deficit-free spacetimes. The 3 
ingoing families were simulated using a finest numerical grid of size 4096 points, with a Courant factor of 
0.1; thus 40960 time steps are required per light-crossing time (for some of the critical solutions presented 
in the next section a grid size of 8192 points was used with a Courant factor of 0.2) For the time-symmetric 
cos^ function we do not need as many points to get good convergence results (because of the milder field 

^though we did check (for the gaussian) that we get the same critical solution when tuning the width, keeping the amplitude 
fixed 
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gradients), so that the highest resolution required for that family was a 1024-point grid. In fact, we get 
acceptable results even after 50 light-crossing times with 1024 points for the cos^ data, whereas the more 
compact ingoing families start having noticeable errors (estimated from convergence tests) in near-critical 
evolution after 3-4 LCT's with 4096 points. 

Fig. 1^ shows the initial scalar field gradient, <i>(r, 0), of typical amplitude for each of the families. Fig. 
H shows the metric function A(r, 0) for a gaussian (the other families have similar shapes for A), and for 
later reference we show how A{r, t) and B{r, t) have evolved at t = 0.6. In order to provide the reader with 
some feeling for the dynamics of a "typical" evolution, Fig. || shows a "space-time" plot of the evolution of 
a sample gaussian with P = 0.1302 that does not form a black hole within 4 LCT's (and it should not, as 
the asymptotic mass of the spacetime is —1.062x10"^). 

3.1 Parameter space survey, varying P 

Figs. ^ and |^ show plots of the asymptotic mass, AI{P), of the spacetime, as a function of the amplitude P, 
about the region M = of parameter space, for the guassian and harmonic families. The second curve on 
each plot shows the initial mass estimate of a black hole (if one formed during the 2 LCTs of the gaussian 
evolution, or 50 LCTs of the harmonic evolution) at the time an apparent horizon is first detected. For 
these amplitudes. Figs. |6| and show the time t and coordinate position r of apparent horizon formation. 
Qualitatively, the features of corresponding plots for the kink and squared guassian (also evolved for 2 LCT's) 
are very similar to those for the gaussian, so for brevity we do not show them. To within the resolution of 
our simulation, the final black hole mass always approaches the asymptotic mass — in other words, we do not 
detect any remnant scalar field (black hole 'hair'). See Fig. || for typical examples. Due to the "reflecting" 
boundary conditions at time-like X, this is not too surprising, although one might have expected something 
like a low amplitude, long wave-length, periodic scalar remnant. The scalar field also tends to zero at late 
times along the event horizon, though in that region of the spacetime our results are not good enough to 
obtain useful decay exponents. 

Fig. 0for the harmonic family shows almost chaotic dependence of the time of AH formation as a function 
of amplitude, as M(P) decreases towards Af = 0. There is evidence that this behavior is also present for 
the other families of initial data, but we have not run those simulations at the necessary resolution to give 
convincing evidence. What appears to be happening is the following. First of all, it is more "difficult" for a 
distribution of the scalar field corresponding to AI > to form a black hole — the distribution needs to be 
compact and centrally condensed. Thus, when we implode a relatively "space-filling" distribution with M 
small (and positive) a black hole will not form on the first bounce through the origin. However, because of 
the boundary conditions at infinity, the scalar field will reflect off T, and, as the field has evolved through a 
strong field (non-linear) regime in the interior, the distribution of energy will be different on the subsequent 
implosion. Moreover, because of the strong gravitational field, the scalar field has a tendency to spend more 
time in the vicinity of the origin on average, preventing it from dispersing throughout the spacetime . So, one 
may expect that if the asymptotic mass M is positive, a region of phase space will eventually be traversed 
during evolution, where it is favorable for a black hole to form, no matter how near-zero is M . However, due 
to the chaotic nature of the curve in Fig. 0, we cannot extrapolate to^M) to io = oo in order to directly test 
this conjecture. 

3.2 The critical regime 

To search for critical behavior in the gravitational collapse of the four families of initial data introduced in 
the previous section, we vary the amplitude P in each case to find the threshold of black hole formation. 
Ideally, we would simply seek the amplitude P* where a black hole forms for P > P*, while for P < P* 
the scalar field bounces around forever without collapse. Unfortunately, such a search is not practical; as 
mentioned in the previous section we do not have the computational resources to follow compact initial data 
for numerous LCT's, and, even with the cos^ data, we do not see any trends that would allow us to conclude 
that if a black hole has not formed after, say, n LCT's, then it probably will not form at all. Thus, what 
we do instead is tune to the threshold of black hole formation on the initial implosion; i.e. we base our 
search on whether or not a black hole forms before any initially out-going radiation reflects off X and then 
falls in, contributing to the collapse. This point of parameter space is labeled as P* in Figs. ^ and ||, and 
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Figure 1: $(r, 0) = (/>,r('', 0) for each family of initial data studied. The three compact families are initially 
ingoing, thus n(r, 0) = $(r, 0), while the harmonic function is time-symmetric with n(r, 0) = = 2/7r, so 
X is at r = 1). 
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Figure 2: A{r, 0) (left-most figure) for a gaussian with P = 0.133051, as obtained by solving the Hamiltonian 
constraint with B{r,0) — 0. This amplitude is used as an example in section 3.3 when we discuss the 
singularity structure, so for reference we also show yl.(r, 0.6) and i?(r, 0.6). Notice in particular how large 
and negative B is towards the origin, indicating that in this region of the grid we are looking at very small 
scales in the problem (the proper circumference element is f = £ta,n{r/£)e^). 



coincides with the place where the initial mass estimate dips to near zero (though for the harmonic data — as 
mentioned in the caption of Fig. ^ — for amplitudes a little larger than P* an apparent horizon first forms 
further out, engulfing the one that is about to form at the smaller radius; see also Fig. 

Near this threshold, it turns out that shortly after the initial implosion, the scalar field and geometry 
close to the origin evolve towards a universal, continuously self-similar (CSS) form. We remind the reader 
that a function which is CSS depends only on a single scale-invariant variable x. Now, the coordinates 
(r, t) in which we solve the equations of motion are not well-adapted to self-similarity. However, after some 
experimentation we found that a natural scale-invariant independent variable in our system is 

x = f, (41) 

where f = ^tan(r/^)e^ is proportional to the proper circumference of an r = constant ring, and tc is 
proper time as measured by the central (r = 0) observer. By convention, tc is negative and increases to 
the accumulation point t* = 0. To better visualize the CSS behavior, we also transform to logarithmic 
coordinates: 

Z = ln(f), T=-\n{tc). (42) 

A CSS function, f{x) = /(e^+"^), then looks like a wave propagating to the left with unit velocity as T 
increases to oo. 

Figs. p|-p^ show scale-invariant functions (j)^z{Z,T), (j)^zz{Z,T) (the second derivative better demon- 
strates the "wave nature" of the critical solution), the mass aspect M{Z,T), and {f'^R){Z,T), for a gaussian 
evolution with P — 0.133059219, which is 'close' to the critical solution (ln(P — P*) = —17.5; see sec. 



3.2.1). In principle, the closer to criticality we tune the initial pulse, the longer the scale-invariant behavior 
should persist in logarithmic space. In practice, of course, finite computational precision and grid resolution 
prohibits fine-tuning to arbitrary accuracy — the figures plotted here show data which is about as close to 
criticality as we can get with 8192 grid points. In terms of the mass aspect in Fig. one can surmise 
that the critical solution is (locally) a kink- like transition from the AdS value M = — 1 to a zero mass state; 
though, interestingly enough, the value of the curvature scalar R at the origin diverges like l/t^ as one 
approaches the accumulation point (we will discuss this in more detail below; also, bear in mind that in Fig. 
O we are plotting f^i?, not R itself). This behavior of the mass aspect suggests that the the transition at 



3 RESULTS 



12 



12.9 



t=3.5 



t=2.5 



t=1.5 



t=0.5 



0.0 



0,0 



--8.9 
1 scale (linear) 



Figure 3: A plot of <P{r,t), the spatial gradient of the scalar field, for sample gaussian initial data with 
P = 0.1302. In this case, a black hole is noi formed. This plot clearly demonstrates the nature of the Dirichlet 
boundary conditions on at X (r = 1 in these coordinates). Even though a black hole does not form, back 
reaction is significant here — notice the non-linear interaction between ingoing and outgoing components of 
the field: when the ingoing and outgoing pulses cross, the ingoing component is amplified, while at the same 
time the outgoing component is surpressed. The effect is most apparent on this plot at around t — 3 near 
the outer boundary; and note that the initial outgoing component of the field is quite small and not visible 
in the picture. 
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Figure 4: Asymptotic mass as a function of pulse amplitude for an initially ingoing gaussian (|2^) of width 
0.05, centered at r = 0.2 in a cosmology with £ = 2/Tr. For the amplitudes that formed an apparent horizon 
within the simulation time of i = 2, the mass estimate at time of AH formation is also shown (it is not clear 
in the figure but this curve does not touch the asymptotic mass curve). The dashed vertical line, labeled by 
P*, is the critical amplitude — see sec. |3.2| . 
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Figure 5: Asymptotic mass as a function of pulse amplitude for the time-symmetric n = \ harmonic function 
(p9[). For the amplitudes that formed an apparent horizon within the simulation time of i = 50, the mass 
estimate at time of AH formation is also shown. The dashed vertical line, labeled by P* , is the critical 
amplitude as discussed in sec. 



3.2 



Notice the discontinuity of the initial mass estimate curve just to the 
right of P* ^ and compare the gaussian case in Fig. ^ The reason for the sudden jump, and difference from 
the gaussian case, is that around t — \ for those amplitudes near P* an apparent horizon is close to forming 
in two locations; to the left of the discontinuity it first forms at larger radii, to the right at smaller (see Fig. 
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Figure 6: The initial coordinate position (r) and time (t) of AH formation for the same set of ampHtudes 
as in Fig. ^ for the gaussian family (if an AH formed within t = 2). 
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Figure 7: The initial coordinate position (r) and time (t) of AH formation for the same set of ai 
as in Fig. ^ for the harmonic familty (if an AH formed within t — 50). 
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Figure 8: Black hole mass estimates as a function of time for gaussians with P = 0,133051 (left) and 
P = 0.1340 (right). The horizontal dashed lines denote the asymptotic masses of the spacetimes. For the 
less massive pulse on the left, the apparent horizon forms after the initial implosion when the field is mostly 
outgoing, and the energy gradually accretes onto the black hole. The more massive pulse on the right forms 
a black hole on the initial implosion, capturing almost all of the scalar field energy except for a small piece 
that initially traveled outwards from t — 0. This piece eventually bounces off T then falls into the black hole. 
Note that the smaller amplitude is less than the critical value, P*, as defined in sec, 3,2, while the larger 
amplitude is super-critical. 



the critical point is Type II — in other words, there is no lower, positive bound on the mass of black holes 
that can be formed by the scalar field. 

Fig, |l^ demonstrates the universality of the solution in the critical regime. Here we plot (j>,zz (as in 
Fig, |l0| for the gaussian) at the same time T for each family in a near-critical evolution. The harmonic 
function appears to have a slightly larger amplitude, but, as we shall now argue, this is apparently just a 



slicing effect. As mentioned in sec, 2,2, because of the gauge that we use, and since we choose to solve for 
A{r, 0) and A^t{r, 0) using the constraint equations (^) and (||), the only slicing freedom we have remaining is 
in the initial conditions for B{r, 0) and B^t{r, 0), Once B{r, 0) and B^t{r, 0) are specified, we have no control 
over the manner in which the slice evolves. For the three compact, ingoing families, the critical behavior 
develops at times ranging from t ~ 0,25 to t = 0,30, and because of the similar initial spatial distribution of 
the energy densities, the evolution has proceeded along very similar slices. On the other hand, the harmonic 
data approaches the critical solution at about t = 1,25, at which point the slice has evolved quite differently 
from the other three families near their respective critical times (we note, however, that by plotting as a 
function of T we do "match" the slices at the origin). To demonstrate that the slices evolve differently, we 
plot in Fig, ^ the normalized inner product between d/dtc (in an f,tc coordinate basis) and Vic for the 4 
families, at the same time used in Fig, |l^. This inner product is the Lorentz gamma factor, W (assuming 
the vectors are time-like), between f = constant observers, and those moving normal to the hypersurface 
tc = constant: 

|a/9ie||vt,| ■ ^ ^' 

This quantity will be the same along identical slices of a spacetime (since such slices will have the same normal 
vectors); thus the harmonic solution slice is clearly different as one moves away from the origin. Another 
interesting feature of this plot for the harmonic data is that it shows gravitational collapse occurring a 
short distance away from the unfolding critical behavior, since, to the right of the peak, the vector d/dtc 
has become space-like (equivalently the surface f = constant has become space-like — see the discussion on 
the singularity structure in sec, [3,3| , and in particular Fig, pO| ), At this point in parameter space for the 
harmonic function there is a lot more mass in the spacetime than that involved in the critical evolution, and 
this is causing an apparent horizon to form at a larger radius (see Figs, || and|^). Presumably, for smaller 
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Figure 9: (j)^ziZ,T) (see ( p2| ) for the definition of Z and T coordinates), for gaussian initial ( p7| ) data with 
P = 0.133059219, a — 0.05, n — 1 and tq = 0.2 in an ^ = 2/7r cosmology. This function of ip is scale- 
invariant in the critical regime, which unfolds roughly between T « 8 and T « 19 (though, interestingly, the 
scale-invariance seems to persist for longer in the scalar field than the geometric quantities — see Figs. |ll| 
and |l|). 
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Figure 10: (j)^zz{Z,T), i.e. the derivative of the function plotted in Fig. ||. 
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Figure 11: The mass aspect, M{Z,T), for the same solution shown in Fig. That M becomes shghtly 
non-monotonic at late times is probably due to numerical error — this is a super-critical evolution, and the 
metric variables are already growing rapidly around T = 19. 
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Figure 12: The Ricci scalar multiplied by (a scale- invariant combination in the critical regime) for the 
same solution shown in Fig. 0. 
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Figure 13: A composite of the scale-invariant function (j>^zz{Z,T) for the near critical solutions of the four 
families of initial data considered, demonstrating universality of the solution in the critical limit. The data 
is plotted at T = 13 (compare Fig. |l^. The harmonic profile appears somewhat different than the others 
due to a slicing effect, as explained in the text (also see Figs. |lj and |l^). See Fig. ^ for the values of P* 
for each family. 



amplitudes one could tune to a threshold solution after several light-crossing times, and perhaps then one 
would more cleanly uncover the critical solution. 

To give more evidence that all the solutions are indeed approaching a universal one in the critical regime, 



we need to compare them on a common spacetime slice. In Fig. 15 we show the same function of the scalar 
field as in Fig. transformed to a Christodoulou type coordinate system (f, w), where a.v — constant curve 
is an ingoing null geodesic We normalized v so that dv = die at the origin; i.e. v also measures central 
proper time. Thus comparing solutions on the same v = constant surface removes any slicing ambiguity |^. 
As can be seen from the figure, the transformed solutions arc all quite similar, though wc lose some accuracy 
in the transformation (which is why we have elected not to use these coordinates in all of the plots in Figs. 

3.2.1 The scaling exponent 7 

Another characteristic feature of Type II critical behavior in gravitational collapse is the universal scaling 
exponent 7 in the relation M = K{p — p*)^'^ . To measure this relationship in the current context, one needs 
to wait for the system to settle down to a steady-state to ensure that the apparent horizon is coincident 



with the event horizon, and hence that the mass estimate (38) gives the correct mass. In AdS, the boundary 
conditions at I prevent us from performing this measurement — initially outgoing radiation that did not 
contribute to the near-critical black hole formation will eventually reflect off X and pollute our measurement. 



However, as discussed by Garfinkle and Duncan |12 , in the near-critical regime (above or below p*) any 



quantity with dimension L', where L is a length scale, should exhibit a scaling relation with an exponent 
of qj. Thus, following those authors, we find the maximum value attained by the Ricci scalar i? at r = 
in sub-critical evolution for tc < 0. Plots of maxt<t^ In |i?(0,t)| vs. ln(P* — P) for the four families studied 

^we are grateful to David Garfinkle for suggesting this procedure to us 
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Figure 14: Tlie Lorentz gamma factor W (|4^ ) between f = constant observer world-lines and those travelling 
normal to the hypersurface tc = constant, at T = 13 for the 4 near-critical solutions as in Fig. |l^. The 
difference between the three initially ingoing families and the harmonic one indicates that we are looking 
at differing slices of spacetime as we move away from r — 0. In fact, the discontinuous peak in 7 for the 
harmonic solution at around Z = —10.2 shows that the f ~ constant surface has become space-like to the 
right of this point, indicating a region undergoing gravitational collapse. 
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Figure 15: The function (j)^zz{Z,v) for the same solutions shown in Fig. v = constant is an ingoing 
null geodesic, chosen here to intersect the origin at T = — Intc = 13 in all cases, and along this hypersurface 
we plot as a function of Z = In f . This coordinate system completely fixes the spacetime slice along which 
we are comparing solutions (at the expense of some loss of accuracy during the transformation), and gives 
additional evidence that there is a universal critical solution. 
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is shown in Fig. Since R oc L~'^, these figures show that the scahng exponent 7 of the 2+ ID AdS 

Klein-Gordon system is about 1.2 ± 0.05. 

Notice that the mass aspect M as defined in (|3^) is dimensionless (which is consistent with the scale- 
invariance of M as plotted in Fig. On the other hand, when we keep i fixed and vary P, the resulting 
black hole mass (being proportional to r^^^) has a length scale of 2, so one would expect the mass-parameter 
scaling relationship for BTZ black holes to go like M — K{P — P*)^"*", where 7 is the same value 1.15 — 1.25 
found above for the scaling of R. The initial-mass estimate curves as shown in Figs. ^ and ^ do roughly 
exhibit this scaling behavior for P > P* . 



3.2.2 Critical behavior in the presence of a point particle 

Here we briefly show how the presence of a point particle (angle deficit) alters the critical solution. The 
particle contributes to the mass of the spacetime (p6[), so the more massive the particle (up to the maximum 
Mpp = 1 in our units) the less scalar field energy is needed to form a black hole, and consequently we have 
smaller amplitudes, P*, at threshold. Interestingly, we find the same critical solution in all cases (see Fig. 
p7| for 3 examples), the only noticeable differences being a systematic phase shift in T related to the mass of 
the particle. The kink-like transition in the mass aspect has the same shape as well, but it ranges from the 
particle mass at r = to M = 0. To within the resolution of our simulations (which was at 2048 grid-points 
in this case) the critical exponent is also the same, namely within the range 7 = 1.15 to 1.25. 



3.2.3 The critical solution from a CSS ansatz? 

Given that we have self-similar behavior in the critical regime, it would be useful to find the exact solution 
assuming a CSS ansatz. Traditionally this is done by assuming the existence of a homothctic Killing vector. 
^ (see 1) 

C^gab = 2gab. (44) 

This implies that in coordinates adapted to the homotheticity, so that ^ = d/dr, each component of gab has 
the form e^"^/, for some function / independent of r. Furthermore, C^Rab = 0, so that Rab and hence the 
Einstein tensor Gab are independent of r. This ansatz is not consistent with the field equations (|^) in the 
presence of the cosmological constant if we assume that the scalar field is self-similar (see [^), as we observe 
in the collapse simulations. Essentially, the scalar field stress-energy tensor (^ would need to decouple into 
a piece that exactly cancels the cosmological constant term plus a scale-invariant term, but we do not think 
that this is possible for a minimally-coupled scalar field. 

It may be that in the 2-l-lD AdS system a different symmetry, such as a conformal Killing vector, would 
be needed to generate the critical solution. Or perhaps the critical solution is only approximately homothetic 
over a limited region of the spacetime. Nevertheless, we have not yet found a symmetry-reduced system that 
reproduces the observed critical behavior.^ 



3.3 Singularity structure 

In all of the solutions that we have studied so far we find that after an apparent horizon forms what appears 
to be a spacelike curvature singularity develops within the horizon. Specifically, the surface of excision 
along which the metric variables A and B and, consequently, the curvature invariants begin to diverge, 
is spacelike. By itself, demonstrating a spacelike surface of arbitrarily large curvature is not sufficient to 
prove that the singularity is spacelike — a counter-example would be the mass-inflation null singularity 
^. However, if we extrapolate to the surface of infinite curvature, based upon the growth of the Ricci scalar 
prior to excision, we still find a spacelike surface (in fact, R grows so rapidly prior to excision — roughly like 
along an r = constant surface if we translate t to zero at the singularity — that the surface of infinite 
curvature essentially coincides with the excision surface at the resolution of Fig |l^ below). In addition, 

^Note added in preparation: David Garfinkle has very recently found a CSS solution in the limit athere the cosmological 
constant vanishes that appears to quite accurately describe the critical solution that we have found p]|]. His result is quite 
intriguing — the cosmological constant is essential for black holes to form, yet apparently it plays very little role in the solution 
at the threshold of formation! 

*We are grateful to Lior Burko for pointing this out to us 
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Figure 16: The maximum of In |-R(0, t <tc)\ as a function of ln(P* — P) for sub-critical {P < P*) evolutions 
of the 4 families of initial data considered. These plots indicate that the maximum of -R(0, t < tc) attained 
during evolution is proportional to (P* — P)~'^^, with 7 Ri 1.15 — 1.25. 
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Figure 17: A composite of the scale-invariant function cj).zz{Z, T) a.tT — 13 for the near critical solutions of 
the gaussian family (cr — 0.05, tq ~ 0.2) with 4 different initial values for A{0, 0) — 0, 1, 2 and 3, corresponding 
to the presence of point particles at the origin with masses 0, 0.864665, 0.981684 and 0.997521 respectively 
(p6|). It is striking that these solutions only differ by a phase in T related to the particle masses; they 
have evolved to the same amplitude after starting with quite different initial amplitudes (namely, P « 
0.13, 0.078, 0.034, 0.013 from A = to 3). 
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B{t,r) — oo along this surface, indicating that the proper circumference measure £ta.n{r/i)e^ goes to 
zero there (see Fig. ^ below). Thus, as with vacuum BTZ black holes, this singularity is crushing^: any 
extended object reaching the singularity is forced to zero proper circumference, regardless of any angular 
momentum or internal pressures that the object might have. 

Figs. |l^ and[l9| are spacetime plots (essentially Penrose diagrams) of $(r, t) and the Ricci scalar R{r,t), 
respectively, for a gaussian initial pulse with P = 0.133051. On the pictures we have superimposed the region 
of trapped surfaces and the inferred event horizon of the space time, found by tracing a null ray backwards 
in time from the place where the AH meets X on the coordinate grid. Fig. ^ show contours of proper 
circumference for the same solution . Th e point P = 0.133051 in parameter space is slightly sub-critical (as 
we have defined criticality, see Sec. |3.2| ) — a black hole forms because the bit of outgoing energy present at 
t — bounces off T and falls back onto the nearly collapsed scalar field, pushing it over the limit. This gives 
us a very clear view of the interior structure; for a more massive pulse the singularity forms shortly after the 
initial implosion, resulting in a thin sliver of an interior in (r, t) coordinates. 

From Fig. |l^ one can see a striking peak that forms in R after the scalar field has bounced through the 
origin and is travelling outwards. In this particular case R has a value of order —10^" in the interior, it then 
grows to order +10* over a very short distance before decreasing to the AdS value of « —15. This 

near-discontinuous behavior in R is characteristic of sub-critical evolutions, and becomes more extreme as 
one nears the critical solution. 

As one approaches the excised space- like surface in Fig. |l^, R starts to grow very rapidly, reaching values 
up to 1 10^* I before excision (this may not be clear on the figure — we chose the gray scale to highlight the 
near-discontinuous behavior in i?). R actually oscillates between large positive and negative values along 
this surface, but our calculations are not sufficiently accurate to conclude that the oscillation is genuine. In 
particular, R is extremely sensitive to the difference 11^ — $^ (see (11)), and 11^ is usually around the same 
order of magnitude as there. We also note that the maximum value attained by R along the excised 
surface becomes smaller towards T. This is to be expected, since in the 2-l-lD system, some scalar field is 
necessary to produce a value of R differing from the AdS value (again, see (|l^)), and as we move towards 1 
along the excised surface there is progressively less scalar field energy remaining. 



4 Concluding remarks 

We have studied black hole formation from the collapse of a minimally-coupled massless scalar field in 2+1 
dimensional AdS spacetime. Outside of the event horizon the spacetime settles down to a BTZ form; in 
the interior a central, spacelike curvature singularity develops. At the threshold of black hole formation we 
find that the scalar field and spacetime geometry evolve towards a universal, continuously self-similar form. 
When a point particle is present at the origin the critical solution is shifted in central proper time by an 
amount related to the mass of the particle. 

By examining the behavior of the curvature scalar during sub-critical evolution we deduced that the 
universal scaling exponent 7 for this system is roughly 1.2 ± 0.05. This value is quite different from the 
scaling exponent 1/2 derived by Peleg and Steif |^ for the collapse of thin rings of dust and by Birmingham 
and Sen [|lOj for particle collisions. However, those works considered different forms of matter, and the phase 
transition was between black hole and naked singularity formation. Thus one would not expect the same 
exponent. Also, the local spacetime geometry about a dust ring or point-particles is necessarily (empty) 
AdS, hence such systems cannot exhibit any of the features, other than mass scaling, that are characteristic 
of critical gravitational collapse. 

Some questions remain unanswered in this work. First, what is the exact nature of the critical solution? 
In other words, what is the character of the symmetry (if any) responsible for the self-similar behaviour, as 
the system does not appear to admit a global homothetic Killing vector |^ Second, will any distribution 

^or deformationally strong, see [|23| . It is straight-forward (though tedious) to see that r = in the non-rotating BTZ black 
hole is a strong singularity as defined by Tipler (though it is not a curvature singularity!). We have not repeated the formal 
calculations in terms of Jacobi fields in our collapse simulations, but because of the central, space-like nature of the singularity 
back-reaction is not likely to weakening it. Note added in revision: shortly after this paper was first published, Lior Burko 
studied the structure of the singularity in 2+lD AdS spacetime using a 'qausi-homogenious' approximation, and did find the 
singularity to be strong and spacelike [ p^ . 



though, as mentioned in the footnote of sec. 5.2.3, David Garfinkle has found a CSS solution that is apparently relevant to 
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Figure 18: The gradient of the scalar field ^{r,t) on the entire solution domain for a gaussian with 
P = 0.133051. On this picture we have also outlined the region of spacetime containing trapped surfaces, 
and drawn in the event horizon with a dashed line (found by tracing a null ray backwards in time from 
the place where the AH reaches r = 1 — which is also T, so our coordinate system breaks down there). We 
stop the simulation at points where the metric variables begin to diverge (the lower boundary of the excised 
region), which presumably is just before a spacetime singularity forms (see Fig ^ for a similar plot of the 
curvature scalar). 
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Ricci scalar 




r 



Figure 19: A plot of the Ricci scalar i?(r, for the same solution as shown in Fig. |T^. During most of 
the evolution \R\ is bounded above by « 10'^'^, but shortly before reaching the excision boundary R starts 
to diverge rapidly, signaling the formation of a spacelike singularity. 
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Contours of constant proper circumference 

1.25 




Figure 20: A contour plot of proper circumference (divided by 27r) f = £tan{r/£)e^ for the same solution 
as shown in Fig. |l^ (the thickness of each contour line is constant in units of proper circumference). This 
plot demonstrates the central nature of the singularity. Along the excised surface approaching T f also 
tends towards zero, though it is not clear with the limited resolution of this figure there. The event horizon 
asymptotes to f = 0.037, i.e. just outside the f = 0.03 contour. 
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of energy that could conceivably form a black hole (i.e. with asymptotic mass M > 0) eventually do so if 
one waits long enough (because of the Dirichlet boundary conditions imposed on the scalar field at X)? A 
third question, related to the first two, is whether the critical solution we have found is a true black-hole- 
formation threshold solution. In other words, that we have a found a universal, CSS solution via a fine-tuning 
process indicates that this critical solution is one-mode unstable; so, does perturbing the critical solution 
"one way" result in a black hole, and perturbing it the "other way" cause the scalar field to remain regular, 
never forming a black hole? The asymptotic nature of AdS spacetime, which is ultimately responsible for 
the boundary conditions of the scalar field at X, prevent us from answering this question in our collapse 
simulations. 

With regards to future work, it would be useful to extend these results to different scalar- field/geometry 
couplings, include a mass and potential terms in the Lagrangian, and to add angular momentum to the 
initial data to study the formation of rotating black holes. It would also be interesting to understand the 
critical behavior in light of the AdS/CFT correspondence. Even though our calculation is purely classical, 
there should be a regime where the classical evolution is a good approximation to the full bulk theory, and 
consequently there should be a dual CFT description of the critical phenomena. 
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